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Abstract 



We consider the long time behavior of the scmidiscrete scheme for the Perona-Mahk 
equation in dimension one. We prove that approximated solutions converge, in a slow 
time scale, to solutions of a limit problem. This limit problem evolves piecewise constant 
functions by moving their plateaus in the vertical direction according to a system of 
ordinary differential equations. 

Our convergence result is global-in-time, and this forces us to face the collision of 
plateaus when the system singularizes. 

The proof is based on energy estimates and gradient-flow techniques, according to 
the general idea that "the limit of the gradient-flows is the gradient-flow of the limit 
functional" . Our main innovations are a uniform Holder estimate up to the flrst collision 
time included, a well preparation result with a careful analysis of what happens at 
discrete level during collisions, and renormalizing the functionals after each collision in 
order to have a nontrivial Gamma-limit for all times. 
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1 Introduction 



The one dimensional Perona-Malik equation is the partial differential equation 

=7rT^^- (x,i)e(o,i)x(o,+oo), (1.1) 

which is usually coupled with Neumann boundary conditions 

u^{Q,t)^u^{l,t)^Q Vi>0, (1.2) 

and an initial condition 

u{x, 0) = uo{x) Vx e (0, 1). (1.3) 
Problem (1.1), (1-2), (1.3) is the formal gradient-flow of the functional 

1 

PM{u) := - / log (1 + ul) dx. (1.4) 
2 io 

The convex-concave behavior of the integrand in (1.4) makes (1.1) a parabolic equa- 
tion of forward-backward type, with forward (or subcritical) regime in the region where 
\ux\ < 1, and backward (or supercritical) regime in the region where Im^^I > 1. 

The analogous problem in two space dimensions was introduced by P. Pcrona and 
J. Malik [18] in the context of image denoising. The rough idea is that small distur- 
bances, corresponding to small values of the gradient, are expected to be smoothed out 
by the diffusion in forward regions. On the contrary, sharp edges should be enhanced 
by the backward nature of the equation in regions where the gradient is large. 

This intuition has actually been confirmed by numerical experiments. There are 
some well known shortcomings, such as the staircasing effect observed in supercritical 
regions, but nevertheless the method reveals some stability, and in any case much more 
stability than expected from a backward diffusion process. 

Equation (1.1) is the prototype of all forward-backward parabolic equations such as 
Ut = {Lp'{urc))x, where (/? is a nonconvex integrand. It is also strongly related to forward- 
backward parabolic equations of the form Ut = {4'{u))xx, where is a nonmonotone 
response function (indeed this is the equation solved by the derivative Ux of solutions u 
of ( 1 . 1 ) ) . Such equations attracted a considerable attention in the last years because they 
are involved in several models, from phase transitions to population dynamic (see [19] 
and the references quoted therein). 

A natural approach to an ill-posed problem is to approximate it by more stable ones. 
Following this idea, several authors proved well posedness results for approximations of 
(1.1) obtained via space discretization [12] or convolution [9], time delay [2], fractional 
derivatives [16, 17], fourth order regularization [3], simplified nonlinearities [4]. A sat- 
isfactory understanding of what happens as the suitable parameter vanishes still seems 
to be out of reach. 

Several papers reported numerical experiments on the Perona-Malik equation in di- 
mension one or two. We refer in particular to [2, 6, 10, 11, 16, 17]. All these experiments. 
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although obtained through different approximation methods, seem to reveal some com- 
mon qualitative features. In particular, the evolution seems to happen in three different 
times scales. We call them "fast time" , "standard time" , and "slow time" , according to 
the terminology introduced in [3] . 

• Fast time. In a time interval of order o(l), solutions of approximated problems 
tend to develop microstructures in the concave region, with fast oscillations be- 
tween very small and very large values of the derivative. This is the staircasing 
effect, which causes an instantaneous drastic reduction of the energy in the back- 
ward regime. From the variational point of view this is hardly surprising, due to 
the concavity of the integrand in that region. More surprising is that this effect 
does not extend immediately to the forward regime, as it could be expected after 
remarking that the relaxation of (1.4) is trivially zero. 

Up to our knowledge, there is no rigorous treatment of this phenomenon. On the 
other hand, the existence of dense classes of smooth solutions of (1.1) (see [13, 14]) 
suggests that it is not reasonable to expect the staircasing effect for all initial data 
with both subcritical and supercritical regions, but at most for "generic" such data. 
This remains a challenging open problem. 

• Standard time. In a time interval of order 0(1), solutions of approximated prob- 
lems evolve in order to reduce the energy in the convex region. Rigorous results 
in this time scale arc known only for the semidiscrctc scheme in dimension one 
(in this paper semidiscrete means discrete with respect to space, and continuous 
with respect to time). In [12] a compactness result was proven, according to which 
solutions of approximated problems converge to something (as the size of the grid 
goes to zero), and all possible limits are classical solutions of (1.1) in the subcritical 
region of Uq. 

The characterization of such limits in supercritical regions remains an open prob- 
lem, as well as any compactness result for different approximation methods or in 
more space dimensions. 

• Slow time. After the second phase of the evolution, the energy has been reduced 
almost to zero, and the solution is close to a piecewise constant function. This is 
consistent with the intuitive idea that piecewise constant functions are stationary 
points of PM{u). Since there is almost no energy left, the evolution slows down. 

Nevertheless, in a slower time scale the plateaus of this piecewise constant func- 
tion tend to move in the vertical direction, with jump points which remain fixed 
in space. The vertical dynamic is nontrivial because neighboring plateaus can 
collide, and actually do collide in a finite time. After each collision at least one 
discontinuity point disappears, and the evolution proceeds as soon as the solution 
becomes a constant and there is nothing else to evolve. 

The aim of this paper is a rigorous analysis of the slow time for the semidiscrete 
scheme in dimension one (we refer to section 2.3 for precise definitions). The first three 
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steps in this direction were done by G. Bellettini, M. Novaga and M. Paolini in [5]. 
First of all, they identified the right time-scale, which turns out to be of the same order 
of the inverse of the grid size (namely 0{ri) if the grid size is 1/n). Secondly, they 
identified the system of ordinary differential equations describing the evolution of the 
plateau heights in the limit problem. Finally, they proved that the rcscalcd solutions of 
the semidiscrete scheme converge to the limit evolution described by that system in the 
half-open interval [0, Tging), where Tging is the life span of the solution of the system. 

The proof of their convergence result is based on the construction of suitable subsolu- 
tions and supersolutions, suggested by a formal development of approximating solutions. 
This method reveals some drawbacks. First of all, it requires some heavy computations, 
which in [5] arc carried out at the expenses of choosing a simplified form of the nonlin- 
earity, a very special sequence of initial data, and Dirichlet boundary conditions. More 
important, it seems quite hard to extend these arguments beyond Tging, namely when 
the interaction of plateaus makes the dynamic highly nontrivial. 

In this paper we overcome these difficulties, and we prove a global-in-time conver- 
gence result (Theorem 3.1). Of course the limit problem is defined by restarting the 
evolution after each collision according to the same rule applied to the new (smaller) 
set of plateaus. 

Before restarting the evolution, it is however necessary to prove that it can be ex- 
tended up to Tsing (included). This fact is quite intuitive, and indeed it has been implic- 
itly mentioned (but not proved) in [5] , when the authors say that the system singularizes 
due to collisions, and not to more strange phenomena. On the other hand, the possibil- 
ity that more than two neighboring plateaus collide in the same time makes this issue 
nontrivial. We overcome this difficulty by proving a 1/4-Holder estimate up to Tging (see 
Proposition 2.1, and Proposition 4.3 for the corresponding estimate at discrete level). 

Then we pass to our convergence result, inspired by the general principle that "the 
limit of gradient-flows is the gradient-flow of the limit" . Since approximating solutions 
are gradient-flows of rescaled approximations of (1.4), it is reasonable to expect the limit 
of the evolutions to be the gradient-flow of the limit energy (in the sense of Gamma- 
convergence). Unfortunately, if we want the limit energy to be finite, we are forced to 
fix a priori the number of discontinuities, and we are back to the interval [0, Tging). 

The first idea is therefore to renormalize the energy after each collision. If we add 
a constant (depending on the grid size) to each approximated energy, then the approx- 
imated gradient-flows do not change, but the hmit energy can be different. This allows 
to iterate the convergence result after each collision, provided that we arrive up to Tging 
(included) with approximating solutions which remain "well prepared", namely close 
enough in many senses to the continuous limit. 

To this end we develop two main tools. The flrst one is a well preparation result 
(Proposition 4.1). When a discontinuity disappears in the continuous limit, then in 
the corresponding interval of the grid the discrete derivative of approximating solutions 
crosses the critical threshold, switching from the concave region to the convex one, 
and instantaneously the discrete solution becomes a well prepared approximation of 
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the new piecewise constant function with a smaller set of plateaus. The second tool 
(Proposition 4.2) is a convergence result up to the first collision time (included), which 
by itself improves the convergence result of [5]. We prove it by rewriting both the 
approximating problems, and the limit problem, in terms of integral inequalities instead 
of differential equations. This formulation, inspired by the theory of maximal slope 
curves (introduced in [8], see [1] for a modern presentation), happens to be much more 
stable when passing to the limit. 

Our techniques work with general nonlinearities (we only need the convex-concave 
behavior of the integrand) , general sequences of initial data (we do not even assume the 
boundcdncss of the energy), and general boundary conditions (we work with Neumann 
boundary conditions because this is the natural choice in applications, but the same 
arguments apply to Dirichlet or periodic boundary conditions). 

Of course several problems remain open. Apart from the notorious questions con- 
cerning fast time and standard time, it could be interesting to prove similar results 
for the slow time in higher dimension, or again in dimension one but with different 
approximation methods. 

A partial contribution in this direction is due to G. Bellettini and A. Fusco [3]. They 
considered a fourth order regularization of (1.1), corresponding to adding a vanishing 
second order term to (1.4). They identified the time-scale of slow time, they computed 
the Gamma-limit of the rescaled energies, and they conjectured that the limit problem 
is the gradient-fiow of the limit energy. Unfortunately in that case this remains a 
conjecture, since up to now no convergence result (even before collisions) is known. 

The limit conjectured in [3] evolves once again piecewise constant functions, but the 
law of the vertical motion is in their case different (the system of ordinary differential 
equations is similar, but with different exponents). This suggests two remarks. On the 
one hand the existence of a slow time vertical motion is a qualitative feature which is 
intrinsic in the nature of (1.1). On the other hand, what exactly happens in the slow 
time from the quantitative point of view does depend on the approximation method. 

This paper is organized as follows. In Section 2 we introduce the rescaled semidiscrete 
scheme, the variational setting, and the limit problem. We also recall the previous 
results which are needed throughout this paper. In Section 3 we state our main results. 
In Section 4 we present the basic tools of our analysis. In Section 5 we collect all proofs. 

2 Notation and definitions 

2.1 Functional spaces 

Continuous setting The more general ambient space we consider is L^((0, 1)), short- 
ened to when it is clear that we are working in the interval (0, 1). We write ||m||lp((o,i)), 
or simply \\u\\p, to denote the p-norm (p e [0, +oo]) of a function u, and {u, v) to denote 
the scalar product of the functions u and v in the appropriate space. 

Let D C (0, 1) be a finite set, and let k :— \D\. The elements of D divide (0, 1) into 
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{k + 1) subintcrvals. 

We call PCd the space of functions which are constant in each subinterval, with the 
agreement that the constant values in any two neighboring subintervals are different. 
In other words, elements of PCd are piecewise constant functions with exactly k jump 
points located in the discontinuity set D. 

We call PSd the space of functions which are Lipschitz continuous in each subinter- 
val, with Lipschitz constant less than or equal to 1, with the agreement that for each 
d & D the limit as x — > d~ is different from the limit as x ^ d'^. In other words, 
elements of PS^ are piecewise suhcritical functions with exactly k jump points located 
in the discontinuity set D. For every u e PSd, and every d & D, the jump height of u 
in d is defined as 

Jd{u) := lim u{x) — lim u{x). (2.1) 

x—^d+ x—>-d~ 

It is easy to see that PCd ^ PSd Q LF'- Every element of PCd is uniquely deter- 
mined by the heights of its (k + 1) plateaus. This correspondence defines an isometry 
between PCd and an open subset of a Euclidean space of dimension {k + 1). When 
needed, we assume that elements of PCd and PSd are defined in the jump points in 
such a way that they are right-continuous. 



Discrete setting Given a positive integer n, we divide [0, 1] into n intervals of length 
1/n, and we consider the space PC„ of all functions which are constant in each subin- 
terval (in this case constants in neighboring subintervals may be equal). The space 
PCn, when endowed with the L^-norm inherited as a subset of L^, becomes a Euchdean 
space isomorphic to M". Since elements of PC„ are thought as L? functions, it is not 
so essential to define them also in points of the form i/n (with i = 0, 1, . . . , n). In any 
case, when needed we assume that the value in any of these points is the same as in the 
interval on its right (on its left in the case i = n). 

Given u e PC^, the discrete derivatives D^^'^u and D~^/'^u are defined as the incre- 
mental quotients 

7-.±l/n / ^ uix±l/n)-uix) ^ ^, 

D'^^''^u{x) := ^ — Vx e [0, 1], 

±l/n 

with the agreement that u has been extended previously to the whole real line (or at 
least to a neighborhood of [0, 1] of width 1/n) by setting u{x) = u{0) for every a; < 0, 
and u{x) — u{l) for every x > 1. 

Given a finite set D C (0, 1) with k elements, we set 

^__^^^rM^_MV[0.1]. (2.2) 

d&D L ^ 

In other words, D„ is the union of all subintervals which intersect D (when d is of 
the form i/n we take the subinterval on its left). It is easy to see that, when n is large 
enough, Dn is the union of k disjoint intervals, and (0, 1) \ D„ has {k + 1) connected 
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components. Since we are interested in passing to the limit as n ^ +00, we can always 
work under this assumption. 

We call PSD,n the set of all functions u e PC„ such that 

\D^/''u{x)\<l^xe[0,l]\Dr,. (2.3) 

The space PSd^h is obviously the discrete counterpart of PSd- In analogy with (2.1), 
the discrete jump height of a function u e PSD,n in a point d ^ D is, defined as 

«(M)-.(HIzi). 

This is equivalent to say that Jd,n{u) := u{x + 1/n) — u{x), where x is any point of 
the subinterval containing d, or of the subinterval on the left oi d ii d — i/n ior some 
i = 1, . . . , n — 1. We point out that (2.4) makes sense for every d, G [0, 1], and not only 
for d E D. Of course we have that \u{x + 1/n) — u{x)\ > 1/n if and only if a; £ Dn- 

The subcritical incremental quotient of a function u e PSD,n is defined as 

Due to (2.3) we have that < SQn{u) < 1 for every u G PSD,n- The subcritical 
incremental quotient is the discrete counterpart of the Lipschitz constant in the intervals 
between discontinuities. 



2.2 Functionals 

A discrete approximation of (1.4) is obtained by replacing the derivative with discrete 
derivatives. Thus we introduce the functionals PM„ : PC„ — > R defined by 

1 

PM„{u) = - / log (1 + iD^/^'uix)]^) dx W G PC„. (2.6) 

The time rescaling due to the "slow time" leads us to consider also the functionals 
nPMn{u). More generally, for each nonnegative integer k we consider the sequence of 
functionals Gn^ : PC„ — >■ M defined by 

Gl!^\u) := ^ / log (1 + \D^/Mx)f) dx - klogn Vti G PC„. (2.7) 

It is clear that Gn\u) (which sometimes we call k-energy of u) coincides with 
nPM„{u) up to an additive constant, and in particular these functionals have the same 
gradient, hence also the same gradient-flow. 

Computing the gradient is a simple exercise in finite dimension. It turns out that 

WGll^\u) = nVPMM = -nD-^l- 



jjl/n 



U 



1 + iDVr^-uP 



"in G PC„. (2.8) 
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The "second order discrete operator" in the right-hand side of (2.8) needs some 
interpretation in the two extremal subintervals, where it requires to compute values of u 
outside [0, 1]. Once again this is done after extending u by setting u{x) = u{0) for every 
a; < 0, and u{x) = u{l) for every x > 1. This agreement is the discrete counterpart of 
the Neumann boundary condition, which in this sense is now included in the right-hand 
side of (2.8). 

Finally, for every finite set D C (0, 1) with \D\ =: k, and every u G PSd, we set 

G^^\u):^J2^og\Mu)\. (2.9) 

deD 

The following result justifies the notation used for G'^\ and shows that the sequence 
Gn\u) has a less trivial hmit as n — >■ -|-oo with respect to the sequence nPMn{u). A 
proof of this result can be found in [7], or simply deduced from the theory of convex- 
concave integrands developed in [15]. 

Theorem A (Gamma-limit of discrete functionals) Let k he a nonnegative inte- 
ger. Let us extend PM^ and Gn^ by setting them equal to +oo for every u E L'^\ PCn- 
Then we have that (all V-limits are intended with respect to L"^ -metric) 

if u is constant, 



and 



V- lim nPMJu) = ■> , 

n->+oo I -l-oo otherwise, 



—oo if u E PCd for some D C (0, 1) with \D\ < k, 
T-Jim^ Gl^\u) = { G'^^(ti) if u e PCd for some D C (0, 1) with \D\ = k, 

-\-oo otherwise. 



2.3 The semidiscrete scheme 

The semidiscrete scheme for the one dimensional Perona-Malik equation is the gradient- 
flow of (2.6). This leads to the problem 

<it) = -VPMn{Vn{t)) Vt > 0, (2.10) 

Vn{^) = U^n, (2.11) 

where {rton} is a suitable sequence of initial conditions with uq^ £ PCn for every n>l. 
The behavior of Vn{t) as n — >■ -|-oo is the subject of the "standard time" theory. 

In the "slow time" theory we speed up the evolution by considering the sequence 
Unit) '■= Vn{nt). It is very simple to show that this sequence solves the rescaled problems 

u'^{t) = -nVPMn{un{t)) yt > 0, (2.12) 
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M.„(0)=Mon. (2.13) 

Thanks to (2.8), the differential equation (2.12) is equivalent to 

u'^it) ^ -VG^^\un(t)) yt>0. (2.14) 

All these problems admit a unique solution defined for every t > 0. Indeed PCn is 
a finite dimensional vector space, the functionals we consider are of class C°°, and their 
gradient (2.8) is globally Lipschitz continuous. It is worthwhile to notice that formula 
(2.8) makes (2.10) the discrete counterpart of (1.1). 

Throughout this paper we consider Un both as a function Un{t) of the time variable 
with values in L^, and as a function Un{x,t) of {x,t) with real values. 

The following properties of u„ are used several times. The proof can be deduced 
from the corresponding properties of Vn stated in [12]. 

Theorem B (Properties of approximating solutions) Let n be a fixed positive in- 
teger, letuQn £ PCn, and letUn [0, +oo) — >■ he the solution of problem (2.12), (2.13). 
Then the following properties hold true. 

(1) (Regularity) We have that 

Un e C°° ([0, +oo); PCn) C C~ ([0, +oo); L") . 

(2) (Standard gradient-fiow estimate) Let k be any nonnegative integer. Then the 
function t — >■ Gn\un{t)) is nonincreasing, and for every < s <t we have that 

\\Un{t) - Un{s)h < {Ci'\un{s)) - Ci'\un{t))f' \t - s\'/'. 

(3) {LP estimate) The function t — )■ ||Mn(2^) ^)||lp((o,i)) ^■^ nonincreasing for every p > 1 
(including p — oo). 

(4) (Total variation estimate) The function t ^ ||D-'-/"M„(a;, t)||j;,i((o,i)) is nonincreas- 
ing. The same is not necessarily true for p -norms with p > 1. 

(5) (Asymptotic behavior) The function Un{t) tends, as t ^ +oo, to the constant 
function equal to the average of uon- 

(6) (Monotonicity and extinction of supercritical regions) Supercritical regions are 
nonincreasing (as set valued maps), and they disappear after a finite time. 

In other words, if uon G PSD,n for some finite set D C (0, 1), then there exist 
j G N, and a finite sequence of times 

= To<Ti< ... <Tj < Tj+i = +00, 

and a finite sequence of subsets 

D = D{0) D D{1) D ...D D{j) = 
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(with strict inclusions) such that 

Unit) e PSD(i),n yt e [Ti,Ti+i), Vi e {0, l, . . 

The fact that supercritical regions disappear in a finite time (which of course depends 
on n) has probably never been stated in the literature, but it is a simple consequence 
of statement (5). In turn, statement (5) follows from three general facts: the average 
of Un{t) is invariant during the evolution, the limit of a gradient-flow is a steady state 
solution, and (2.8) is zero if and only if u is constant. 



2.4 The limit problem 

Let k he & nonnegative integer, and let D C (0, 1) be a finite set with \D\ — k. Given 
an initial condition uq e PCd, we define an evolution u{t) starting from uq according 
to the following algorithm. 

If /c = 0, then the initial datum uq is constant, and we define u{t) as the stationary 
solution u{t) = Uq. 

If A; > 0, let D = {di, ...,dk} with < di < . . . < < 1, let do 0, dk+i := 1, 
and let (with i — 0, . . . , k) denote the constant value of Uq in the interval {di, di^i). 
Let (ao(t), ai(t), . . . ,ak{t)) be the (unique) solution of the system of {k + 1) ordinary 
differential equations 

1 1 

di-do ai{t)-ao{t)' 



di+i - di \ai+i{t) - ai{t) ai{t) - ai-i{t) J 
1 1 

dk+i — dk a,k{t) — ak-i{t) ' 

with initial conditions aj(0) = a^i for every i = 0, . . . , k. 

Let u{t) be the piecewise constant function whose value in (dj, dj+i) is aj(t), defined 
as soon as the solution of the system exists. To this end, we have the following result 
(the proof is given is section 5.3). 

Proposition 2.1 Let k, D, uq, u{t) be as above. Then we have the following conclu- 
sions. 

(1) (Local but not global existence) The system of ordinary differential equations has 
a local solution defined on a maximal interval [0,Tsing) with Tging € (0, +00). 

(2) (L°° estimate) We have that 

\Ht)\\oo < IImoIIoo yt e [0,T,i„g). (2.15) 



a'oit) = 
a'kit) = 
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(3) (1/4- Holder continuity up to collision) For every {s,t) e [0, Tging)^ we have that 
hit) - uis)\\2 < (Skf exp (J^G^^\uo)^ \t-s\'/\ (2.16) 

Thanks to the Holder continuity (2.16), we can define u{t) up to Tging (included). This 
extension fulfils (2.16) and (2.15) in the closed interval. 

Moreover, u{t) e PCd' for some D' CD. li D' — D we can continue the solution 
of the system of ordinary differential equations beyond Tging, but this contradicts the 
maximality of Tging. Thus D' is strictly contained in D, which means that at time 
t = Tsing we have a collision between at least two adjacent plateaus. At this point we 
restart the construction of u{t) from M(Tsing), which has a smaller set of jump points. 

This procedure defines a function u e C° ([0, +00); L^) with u{0) = uq. For each 
t > the function u{t) is piecewise constant in the space variable, and its discontinuity 
set is contained in the discontinuity set of uq. There is a finite set of singular times when 
two or more adjacent plateaus collide, hence one or more discontinuities disappear. After 
each collision, the involved plateaus remain attached forever, and the evolution goes on 
according to the same rule apphed to the new set of plateaus. After the last collision 
u{t) becomes constant, and it does not move any more. This constant is actually the 
average of uq (indeed the average of u{t) is invariant during the evolution). The function 
u is of class with respect to the time variable outside the finite set of collision times, 
and uniformly continuous as a function from [0, +00) to L^. 



3 Main results 

The main result of this paper is the following convergence result. 

Theorem 3.1 (Global-in-time convergence) Let D' Q D Q [0, 1) be two finite sets. 
Let Uq G PCd', and let {won} be a sequence such that 

Uon e PSD,n Vn > 1, (3.1) 

uon^uo mL2((0,l)). (3.2) 

For every n > 1, let Un(t) be the solution of the approximating problem (2.12), 
(2.13). Let u{t) be the solution of the limit problem with initial condition uq, as defined 
in section 2.4- 

Then we have the following conclusions. 

(1) (Global- in-time L^-convergence) We have that Un{t) u{t) m C° ([0, +oo); L^), 
namely 

lim sup \\un{t) - xi(t)||L2((o,i)) = 0. (3.3) 
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(2) (Global-in-time "uniform" convergence) For every t > 0, let D{t) be the disconti- 
nuity set ofu{t), and let Dnit) he the union of all subintervals containing elements 
of D{t), defined according to (2.2) with D(t) instead of D. Let us set 

Kn := {(x, t) e [0, 1] X [0, +oo) : a; ^ D^{t)}. (3.4) 
Then we have that 

lim \\un{x, t) - u{x, t) \\L-o(Kn) = 0. (3.5) 

n—>-+co 

We considered a piecewise constant initial datum Uq because this is the natural space 
where the Gamma-limit of the renormalized functionals is finite. On the other hand, 
we emphasize that the approximating sequence {won} is quite general. In particular, we 
did not assume that it is a recovery sequence, or that its energy is bounded, and the set 
of discrete jump points of mo„ can be larger than the set of jump points of Uq. What is 
essential is that all discrete jump points of approximating functions are contained in a 
fixed finite set D (if not, there are counterexamples even to local- in-time convergence). 

Another reason for looking at piecewise constant data is that they are expected to 
be the limit as t — > -|-oo of evolutions in "standard time" . Up to our knowledge, this 
has been proved rigorously only for "generic" piecewise subcritical data (see [6, 12]). 

Whenever the "standard time" evolution admits a piecewise constant limit, we can 
start our "slow time" analysis from that limit. We state this idea formally in the next 
result. We point out that in this case we do not assume that the sequence {won} of 
initial data has a piecewise constant or piecewise subcritical limit (actually we do not 
even assume that it has a limit) . 

Theorem 3.2 (Convergence for more general initial data) Let D C (0, 1) be a 

finite set, and let {uon} be a sequence satisfying (3.1). 

For every n > 1, let Vn{t) be the solution of problem (2.10), (2.11) (no time rescal- 
ing), and let Un{t) = Vn{nt) be the solution of the rescaled problem (2.12), (2.13). Let 
us assume that there exist 5* > 0, f G C°([5', +oo), L^), D' C D, and v^o G PCoi such 
that 

lim Vn{t) = v{t) Wt > S, (3.6) 

n— ^+oo 

lim v{t) = Voo, (3.7) 

t—^+oo 

where both limits are intended in L^. Letu{t) he the solution of the limit problem, defined 

as in section 2.4, with initial condition ?i(0) = v^. 

Then for every T > we have that Un{t) u{t) in C° ([T, -|-oo); L^), namely 

lim sup \\un{t) - w(i)||L2((o,i)) = 0. (3.8) 

n^+oo 

Moreover, if fCn is defined as in (3.4), we have that 

lim \\un{x,t) - M(x,i)||ioo(x;„n([o,i]x[r,+oo))) = 0. (3.9) 
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We conclude by mentioning a possible extension of our results. 

Remark 3.3 For the sake of simplicity, we devoted this paper to the model case of 
the Perona-Malik equation, in which the integrand is V'(o") := log(l + o"^). Similar 
techniques apply to larger classes of convex-concave integrands, for example the case 
where v?(cr) := a~^{l + a'^)"/^ for some a E (0, 1). In this case the "slow time" is of 
order 0{n}~°'), the hmit energy is 

Ga{u):=J2\M^T^ (3-10) 

and the system of ordinary differential equations governing the evolution of the plateau 
heights is 

,/ V _ 1 / ai+i{t) - aijt) ai{t) - ai^i{t) \ 

""'^ ' ~ d^+l (k V - a,(t)|'-" \a,{t) - a,_i(t)|'-V ' 

suitably modified for i = and i = k. 

There are, however, some remarkable difi:erences. On the one hand, this situation is 
simpler because the hmit energy (3.10) is bounded from below, hence there is no need 
to renormalize it after each collision. On the other hand, the limit energy can be finite 
even if u has infinitely many jump points. 



4 Fundamental tools 

In this section we state the main ingredients needed in the proof of our main results. 

The first one is a well preparation result, which plays its role at the beginning of 
the evolution and during each collision. In input we have a sequence satisfying (3.1) 
and (3.2) as in the assumptions of Theorem 3.1. In particular, some of the jump points 
might disappear in the limit (this happens if and only if D' is strictly contained in D), 
and there is no information on the /c-energy or the A;'-energy of the sequence (where 
k := \D\ and k' := \D'\). We prove the existence of a sequence of times 5*^ — )■ such 
that Un{Sn) is a "well prepared" sequence, namely it still converges to Uq, all its elements 
lie in the corresponding space PSD',n, ^-nd their /c'-energies converge to the k' -energy of 

Proposition 4.1 (Well prepeiration) Let D, D' , {u^n}, Uq, Un{t) he as in Theo- 
rem 3.1, and let k' :— \D'\. 

Then there exists a sequence Sn ^ of positive real numbers such that 

UniSn) £ PSD',n foT cvcry u large enough, (4.1) 

lim Gf)(«„(5„))=G(^')(«o), (4.2) 

n— ^--l-oo 

lim max \\unit) — Uq\\2 — Q ■ (4.3) 

n->+oo te[0,5„] 
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The second tool is a convergence result up to the first jump extinction. It plays its 
role in the time intervals between collisions. Now in input we have a "well prepared" 
sequence, or at least a sequence of initial data with bounded energy, and without vanish- 
ing jump points. We prove some sort of uniform convergence on an increasing sequence 
of time intervals (depending on n). As n — >■ +00 these intervals invade the whole time 
interval up to the first collision. 

Proposition 4.2 (Convergence up to first jump extinction) Let k be a positive 
integer, and let D C (0, 1) he a finite set with \D\ — k. Let uq G PCd, and let {uon} be 
a sequence satisfying (3.1), (3.2), and 

supG'(f)(iion) < +00. (4.4) 

n>l 

For every n > 1, letUn{t) be the solution of the approximating problem (2.12), (2.13), 
and let 

Tsing.n := sup {t > : M„(t) e PSd,.} (4.5) 

be the first time when a discrete jump disappears (it is the time Ti in statement (6) of 
Theorem B). Let u{t) and T^^g be defined as in section 2.4. 
Then there exists a sequence {Tn} of real numbers such that 

< r„ < Tsing.n Vn > 1, (4.6) 
lim T^^Tsing, (4.7) 

n— >'+oo 

lim max \\uJt)-u(t)\\2^0. (4.8) 

Moreover we have that 

lim G^^\un{t)) = G^^\u{t)) yt e (0,T,i„g). (4.9) 

Finally, we present a qualitative property of approximating solutions which could 
be interesting in itself. It is the discrete analog of statement (3) of Proposition 2.1. 
We point out that the standard gradient-fiow estimates (statement (2) of Theorem B) 
control the 1/2-Holder constant of Un{t) in terms of the descent of the energy, but such 
estimates are useless if the energy is not bounded from below independently on n, and 
this is exactly what happens in this model when t approaches a collision time. 

The following 1/4-Holder estimates overcome this difficulty. 

Proposition 4.3 (Uniform l/4-H61der continuity) Let k, D, {uon}, Un{t), Tsing,n 
be as in Proposition 4-2. 

Then for every n>l and every {s,t) e [0,rsing,n]^ we have that 

\\u„{t) - Ur,{s)h < (3^)3/4 exp (^JLG(f)(«on)) \t-s\'/\ 
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5 Proofs 



5.1 Basic estimates 

In this section we collect some general facts, which are going to be used several times 
in the proof of our main results and basic tools. The first one concerns "double limits" 
(we omit the simple proof). 

Lemma 5.1 (Double index sequence) Let {Am,n} (with {m,n) G N"^) be a double 
indexed sequence with values in a metric space X. Let us assume that for every m G N 
there exists 

■^m.oo ■ — lim Ajjin, 

n—^+oo 

and that there exists 

A-r^ r«n llUl A. 



m.oo- 



Then we have the following conclusions. 

(1) (Standard conclusion) There exists a sequence ruk — >■ +oo of nonnegative integers 
such that 

lim A^^^k = ^oo,oo- (5.1) 

K->+0O 

(2) (Refined conclusion) For every sequence of real numbers r^ — >■ +oo there exists a 
sequence ruk — )> +oo o/ nonnegative integers such that ruk < for every k large 
enough, and such that (5.1) holds true. □ 

In the next result we estimate from below the norm of a discrete derivative. In the 
continuous setting, when we know the values of some function / G -f^o((0' 1)) some 
given points, then we can estimate ||/a;||2 from below. The conclusions of the following 
lemma are the discrete counterpart of such estimates. 

Lemma 5.2 (Discrete derivative estimates) Let n be a positive integer. Let f G 

PCn be a function which is equal to in the last subinterval (1 — 1). Let us consider 
the discrete derivative D~^/'^ f{x), defined after setting f{x) = in (— 1/n, 0). 
Then the following estimates hold true. 

(1) We have that 

\\D-''-f{x)\\^>2\\f{x)\U. (5.2) 

(2) Let k be a positive integer, and let < di < . . . < dk < 1. Then we have that 

with the agreement that do — 0, d^+i = 1, and f{do) — f{dk+i) — 0. 
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Proof Let fi (with i = 1, . . . ,n) denote the value of / in the interval {{i — l)/n,i/n). 
Let us set /o := (this choice is consistent with our extension of f{x) in (— 1/n, 0)), and 
let us recall that fn — due to our assumption on /. Let j e {1, . . . ,n} be the index 
(or one of the indices) such that ||/(x)||oo = \ fj 
Then we have that 



i=j+l 



1=1 



> 



E(/. - f^-l) 



i=l 



+ 



\fj-M + \fn-fj 



2|/,|, 



i=j+i 

which proves (5.2). 

Let us consider now the second statement. Let us set io := 0, ik+i '■= n, and 
ih '■= \ndh\ + 1 for every h = 1, . . . ,k. With this notation we have that f{dh) — fi^^ for 
every /i = 0, 1, . . . , + 1. Moreover it is easy to see that 

ih+i ih \ndh+i\ + 1 \ndh\ + 1 ^ , j , 1 
— S "^+1 — «/i H — ■ 

n n n n n 

Thus from Holder's inequality it follows that 



|D-^/"/(, 



X 



L2((0,l)) 



h=0 
k 

s E 

h=0 
k 



[ih+i/n) - (ih/n) 

> lP"'/"f( 



D ^^"/(a:^)llLi((ift/„,i;,+i/n)) 



'^^i(.ih/n,ih+i/n))- 



(5.4) 



To be precise, the first inequality requires that all indices ih are distinct, and this 
is true only when n is large enough. On the other hand, the final result is true in any 
case, because it is enough to ignore the intervals of length zero in the first sum. 

Finally we have that 



iD-'/^'fix 



L^{{ih/n,ih+i/n)) 



E (/^ - /-i) 

i=ih+i 



E \f^-M> 
i=ih+l 

= l/ih+i - /iJ = l/(4+i) - /(4)|- 

From the last estimate and (5.4) we obtain (5.3). □ 

In the next statement n is fixed, and we present several estimates relating fc-energies, 
the norm of their gradient (the slope), jump heights, and subcritical incremental quo- 
tients. These estimates are the technical core of our analysis. 
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Lemma 5.3 (Fundamental estimates) Let k and n be positive integers, let D C 
(0, 1) be a finite set with \D\ = k, and let v G PSD,n- Let Dn be defined according to 
(2.2), and let us assume that n is big enough so that (0, 1) \ Dn has {k + 1) connected 

components. 

Let Gn\v) be the functional defined in (2.7), let VG'i'^^(f) be its gradient, let Jd,n{v) 
be the discrete jump heights of v defined in (2.4), and let SQniy) be the subcritical 
incremental quotient of v defined in (2.5). 

Then we have that 

Gi'n^)>^log(^mm|J,»|), (5.5) 
Gl!:\v) < ^log (1 + [SQMf] + i ^log (1 + [J,,Mf) , (5.6) 

deD ^ ^ 

\\VG^:Kv)l>n\SQM\, (5.7) 
\\^Gi'\v)\\^>(unn\Ja,n{v)\^ ■ (5.8) 
Finally, if D = {di, . . . , dk} with < di < . . . < dk < 1, then we have that 



with the agreement that do — 0, d^+i — 1, and Jdo,n{v) — Jdk+i,n{v) — 0. 



Proof of estimates on the functional From the definition of jump lieiglits we liave tliat 
n 

I log [1 + \V^'"'v[x)\^) dx - klogn = - 

deD 



I \og{l + \D^lM^)\'') dx-k\ogn = ]^Y.^og{l + n^JdA'^)f)-^\ogn'' 



deD ^ 



+ [J,^n{v)f]. (5.10) 



In order to prove (5.5), we estimate tfie rigfit-fiand side of (2.7) from below by 
considering only the integration over From (5.10) we deduce that 

Gl!:\v)>^ I \og{l + \D''M^)\')dx-k\ogn = \Y.^og(\ + [JaMf) 

JD„ / 

> E^°g l'^<^.'»(^)l - ^^og I min I Jd,n(w) I I , 



deD 

which proves (5.5). 
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On the other hand, from (2.5) we have that 

^ / log (1 + 1^1^(0;) H dx<^log{l + [SQn{v)]'). (5.11) 



'[0,1]\-D„ 

Summing (5.11) and (5.10) we obtain (5.6). 
Proof of estimates on the slope Let us consider the function 



L)i/"t;(a;) 

1 + [Dy^v{x)y 



m-.-n ^ , (5.12) 



It turns out that f{x) = in the last subinterval (1 — 1/n, 1). Moreover VGn\v), 
whose expression has been computed in (2.8), coincides (up to the sign) with the discrete 
derivative D~^^"- f(x), computed after setting f{x) = in the interval (— 1/n, 0). 

Therefore, applying (5.2) to the function f{x) defined in (5.12), we obtain that 

l|VG«(^)||, = \\D-yy{x)\\^ > 2\\fix)\U (5.13) 

Let us estimate the right-hand side of (5.13) from below by restricting the L°°-norm 
to [0, 1] \ D„. Since the function a — > a{l + a'^)~^ is increasing in [—1,1], and since 
< SQn{v) < 1, we obtain that 

II n/ Ml I n\D^/"-v(x)\ ■, X ^ 1 nSQn(v) n , , 

\\f{x)U > max ' \ ^' : X e [0, 1] \ U \ , > - SQM, 

[ 1 + [Dy^'vix)] J 1 + [SQn{v)] 2 

which proves (5.7). 

Now let us estimate the right-hand side of (5.13) from below by restricting the L°°- 
norm to D^- We obtain that 

,1, . , n\D^/Mx)\ \ n^\JdAv)\ 

||/(x)||oo > niax< — - — — 2 : X e -D„ > = max- 



1 + [Dy^v{x)Y J deD 1 + ^2 [J^^^^v)] 



2 ■ 



Since n|Jd^„(f)| > 1 for every d E D, and since the function a — )■ a{l + o"^)~^ is 
decreasing for a > 1, we have that the maximum is achieved when the argument is 
minimum, hence 

n^lJw (v)\ 2n'^ jmn \Jd,n{v)\ 

||VGi'=H^^)||,>2||/(a:)||oo>2max ^l^'^-^^^l - 



deD l + n2[Jd»]' 1 + ^2 



. 2 ■ 

min \Jd,niv) 



Recalling once more that the minimum is greater than or equal to 1/n, estimate 
(5.8) follows. 
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A more refined estimate, keeping into account all jump heights, follows from (5.3) 
applied to the function f{x) defined in (5.12). Since D^^^v{d) = nJd,n{v) for every 
d E D, we obtain exactly (5.9). □ 

In the next lemma we consider the difference between two piecewise constant func- 
tions V and w. The main idea is the following. If the number of discrete jump points is 
finite, and their location is fixed, then the L^-norm ofv — w estimates both the L°°-norm 
oi V — w, and the difference between jump heights. 

Lemma 5.4 (Uniform and jump-height estimates) Let D C (0, 1) be a finite set, 

let Kq be the minimum of the lengths of the (|D| + 1) intervals into which (0, 1) is divided 

by D, and let n > SK^^ be a positive integer. 

Let D' (1 D and D" (1 D be two subsets, and let v G PSD',n md w G PSD",n be two 
piecewise constant functions with discrete jump heights Jd,n{v) and Jd,n{w), respectively 
(we can think both jump heights as defined for every dE D). 

Then we have that 

min [Kq, \\v — < 311^ — w||2^^, (5-14) 

min {Ko, \Jd,n(v) - Jd,nH\} < 6|k - '^WT Vd G D. (5.15) 



Proof Let f{x) := \v{x) — tv{x)\, and let fi (with i = 1, . . . ,n) denote the value of / 
in the i-th subinterval. Let j G {1, . . . , n} be the index (or one of the indices) such that 

ll/(^)l|oo = /.-. 

For every x G [0, 1] \ we have that \D'^/''v{x)\ < 1 and \D^^''w{x)\ < 1, hence 

\D'/^f{x)\<2 yxe[0,l]\Dn. (5.16) 
Let us set H :— mm{Ko, fj}, and let us consider the two intervals 

n 3 ^ n J ^ \ n ' n 3 

Since n > SK^^, estimate (5.16) imphes that in at least one of these intervals the 
difference between the values of / in any two neighboring subintervals is always less 
than or equal to 2/n. Let us assume, without loss of generality, that this happens in / 
(the other case is specular). Then we have that 

fix) > fj-2(^^-x^ > H-2j > I VxG/, 

so that 
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and therefore 

mm{Ko, \\v - w\\oo} ^ H < ?,\\v - ^11^2^) < 3||f - ^y||i2^((o,i))- 
This proves (5.14). Now we have that 

\Jd,n{v) - Jd,n{w)\ = \ Jd,n{v " HI < '^h " ^lloo G D, 



hence 



min {Kq, \ Jd,n{v) - Jd,n{w)\ ] < 2min{i^o, ||^^ - w\\oo} < 6||v - 



w 



i2/3 



for every d e D, which is exactly (5.15). □ 

In the last lemma we consider a sequence Vn ^ v. We point out that this sequence 
is allowed to loose jump points in the limit. 

Lemma 5.5 (Jump convergence and BV estimate) Let D' 'Z D C. (^0,1) be two 
finite sets. Let v e PSd', and let {vn} he a sequence such that Vn ^ v in L^. Let us 
assume that for every n > 1 we have thatVn & PSD"{n),n for some finite set D"{n) C D. 
Then we have the following conclusions. 

(1) (Jump convergence) Let Jd{v) be the jump heights of v, and let Jd,n{'^n) be the 
discrete jump heights of v„ ( we can think both jump heights as defined for every 
d & D, with the agreement that Jd{v) — when d & D \ D'). Then we have that 

lim Jd,nM^Jd{v) WeD. (5.17) 

n— >+oo 

As a consequence, D' C D"{n) for every n large enough. 

(2) (Uniform BV estimates) We have that 

sup {||Di/X||i + ll^nlloo} < +00. (5.18) 
n>l 

Proof of statement (1) Let Wn G PCn be the piecewise constant approximation of v 
defined by 

Wn{x) := V (^^^^ e [0, 1], Vn > 1. 

Let Kq be as in Lemma 5.4, and let n > ?>Kq^ as in that lemma. It is not difficult 
to see that 

\Jd{v)-JdAwn)\<'2n-^ ydeD, 



hence 



\Jd{v) - Jd,nM\ < \Jd{^) - Jd,n{^n)\ + \Jd,n{'^n) - Jd,n{^r, 
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Therefore, applying (5.15) with v = Vn and w = Wn, we obtain that 

min {Kq, \Jd{v) - Jd,n{'"n)\} < min {Kq, \Jd,n{Wn) - Jd,n{'"n)\ + 2n~^} 

< min{Xo, \Jd,n{'^n) - Jd,nM\} + "^n^^ 

< 6\\Wn - Vn\\l^^ + ^n"^ 

< 6 {\\wn - v\\2 + \\v - VuhY^^ + '2n~^ 

for every d & D. All the terms in the right-hand side tend to as n — > +00. This 
proves (5.17). 

Proof of statement (2) Let D„ be defined by (2.2). For every d E D we have that 
D^/'^Vn{x) — nJd,n{vn) fo^ Gvery x in the corresponding subinterval, hence 

ll^'^"^n||i= / \D'/^Vn{x)\dx+ [ \D'/^Vn{x)\dx<J2\Jd,nM\ + l. 

Due to (5.17) we have that 

lim V \Jd,n{Vn)\ + 1 = \ Jd{v)\ + 1 < +00, 

n— >-+oo ' ' 
de-D d&D' 

which is enough to prove the equi-boundedness of the total variations 

The uniform bound on ||f„||oo follows from the uniform bound on total variations, 
and from the fact that the average of Vn tends to the average of v (hence averages are 
equi-bounded) . □ 

5.2 Evolution problems as maximal slope curves 

Both the differential equation (2.12), and the system of ordinary differential equations 
introduced in section 2.4, are equivalent to suitable integral (in)equalities. This equiv- 
alence is the key point in the theory of maximal slope curves, for which we refer to [1]. 

For the sake of simplicity, we want to keep this paper as independent as possible of 
the general abstract theory of maximal slope curves. For this reason, in Proposition 5.7 
we state the two implications we need throughout this paper, and we provide a self- 
contained and almost elementary proof of them. Before stating these implications, we 
need the following definition. 

Definition 5.6 (Slope of limit functional) Let D C (0, 1) be a finite set with \D\ = 
k. Let us assume that D = {di, . . . , dk}, with < di < . . . < dk < 1. Let v € PCd, and 
let Jd{v) be its jump heights. Let be the functional defined in (2.9). 
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The slope of Goo in the point v with respect to the L^-metric is the nonnegative real 
number whose square is given by 



There are several interpretations of the slope. In this case the domain where the 
functional is finite is isometric to an open subset of W''^^, and under this isometry the 
functional can be identified with a function of (k + 1) real variables. The slope coincides 
with the norm of the gradient of this function. 

The approximated solution Un{t) is the gradient-fiow in of the functional Gn^ for 
every t > 0, while the function u{t) defined in section 2.4 is the gradient-fiow in of the 
functional up to Tging. In any case, what we need in this paper (and in particular 
in the proof of Proposition 4.2) are the following two facts. 

Proposition 5.7 (Differential equations vs maximal slope curves) The approx- 
imating problems and the limit problem can be reformulated as follows. 

(1) (Approximating problems) Let n he a positive integer, and let uon G PCn- Let 
Un e ([0, +oo); L^) be the solution of (2.12), (2.13). Then we have that 



for every < s <t and every k eN. 

(2) (Limit problem) Let k be a positive integer, let D C (0, 1) he a finite set with 
\D\ = k, and let uq G PCd- Let Tq > 0, and let v e ((0, Tq); L^) he a function 
such that v{0) — Uq, v{t) e PGd for every t e [0, Tq], and 



for every < s < t < Tq. Then v{t) coincides in [0,7o] with the function u{t) 
defined in section 2.4. 

Proof From equation (2.12), which is the same as (2.14), we have that 





G^^\v{s))-G^^\v{t))>lJ^ \\v'{T)\\ldT+lJ^ \\VG^^\v{T))\\ldT (5.20) 



d_ 
di 



Gl^\u^{t)) = -{VGi^\u^{t)),u'^{t)) 



<mi + l\H'\u^{t))\\i. 



Integrating in [s,t\ we obtain (5.19). 
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Let us consider now inequality (5.20). Let D = {rfi, . . . ,dk} with < rfi < . . . < 
dk < I, let do := and dk+i ■= 1, and let us identify v{t) with the vector of plateau 
heights (ao(t), ai(t), . . . , afe(t)), where ai{t) is the constant value of v{t) for a; e {di, dj+i). 

The if^ regularity of v{t) implies regularity of all components. Let us compute 
the time derivative of the function t G^\v{t)). Using the chain rule, and rearranging 
the terms, for almost every t e (0, Tq) we obtain that 

k 



dt ' dt 

1=1 



^ 1/1 1 \ 

with the agreement to neglect the two fractions involving indices less than or larger 
than k (which appear in the terms of the sum corresponding to i = and i = k). 

Applying the inequality xy < 2~^{x'^ + y^) to each term of the sum, we find that (for 
shortness' sake we drop the dependence on t in the right-hand side of the first line) 

-^G^:^\v{t)) < ^E(rfm-ci.)[a:f + ^E;r^f;r^-;r^)' 

dt 2-^ 2 di+i - di \ai+i - Qi ai-ai^x) 



1 ,2 



2M^v-.„2-^||VG(^)(^(t))||^ (5.21) 
for almost every t e (0, Tq). On the other hand, from (5.20) we easily obtain that 

- I^^H^W) > + \ ||VG^n^W)||^ (5.22) 

for almost every t e (0, Tq). Comparing (5.21) and (5.22) we deduce that there is 
equality for almost every t e (0,To). But in the inequality used to deduce (5.21) we 
have equality if and only if 

1/1 1 



\/ di+i - di a'i{t) 



\ai+i{t) - ai{t) ai{t) - ai-i{t) 



for every i — 0, . . . ,k and for almost every t e (0, Tq), which is equivalent to the system 
of ordinary differential equations introduced in section 2.4. Since the right-hand side is 
continuous, we can conclude that actually ai{t) is of class C^, and we have equality for 
every t in the closed interval [0,To]. □ 

The formulation in terms of differential inequalities is very stable when passing to 
the limit. In the proof of Proposition 4.2 we obtain (5.20) by passing to the limit in 
(5.19) as n — )■ -l-oo. The following result is fundamental in that stage. 
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Proposition 5.8 (Bounded slope sequences) Let k be a nonnegative integer, and 
let D C (0, 1) he a finite set with \D\ = k. Let v G L"^, and let {vn} be a sequence such 
that Vn e PSD,n for every n>l, and such that Vn ^ v in L^. 

Let Jd,nivn) denote the discrete jump heights ofvn, and let us suppose that there exist 
Co > and no > 1 such that 

\Jd,n{^n)\ > Co yde D, Vn > no- (5.23) 

Then we have the following conclusions. 

(1) (Gamma- liminf inequality for slopes) We have that 

]im|nf ||VG(f)(^J||, > pGt\v)\\^, (5.24) 

where the right-hand side is intended to he +00 if v ^ PCd- 

(2) (Bounded slope sequences are recovery sequences) If in addition we assume that 

sup||VG'«(^„)||2<+oo, (5.25) 

n>l 

then we have that v e PCd, and moreover 

\Jd{v)\>CQ WeD, (5.26) 
lim Gi'\vn) = G^^\v). (5.27) 

n— >+oo 

Proof We prove the two statements in reverse order. 

Statement (2) Let Ci denote the supremum in (5.25). Prom (5.7) we have that 

ci > \\VG^^\vn)\\^ > nSQnivn) Vn > 1. (5.28) 

Let 5 > 0, and let Ds denote the neighborhood of D with width 5. Then (5.28) 
implies that 

D^/'^Vnix) uniformly in [0, 1] \ Ds, 

which in turn implies that v{x) is constant in each connected component of [0, 1] \ D^. 

Since d is arbitrary, this proves that v G PCd' for some D' C D. On the other hand, 
assumption (5.23) and the convergence of jump heights (5.17) imply (5.26), which proves 
also that actually D' = D. 

In order to prove (5.27), we split the integral in (2.7) into an integral in and an 
integral in [0, 1] \ For the second one we apply (5.28) and we deduce that 

0<^ / log(l + |i5Vx(a;)|2) da;<^log(l + |5g>„)|2) <^logfl + 4V 
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which proves that the integral in [0, 1] \ Z}„ tends to 0. 

For the integral in D„ we apply (5.10) to the function Vn, and we deduce that 



/ log (1 + \D'/^Vn{x)\^) dx-klogn^ ^ S^log + [JdMf] . 



d&D 

From the jump convergence (5.17) it follows that this expression tends to 

deD 

which completes the proof of (5.27). 

Statement (1) Let us take any subsequence (not relabeled) which realizes the lim inf 
in (5.24). We can assume that (5.25) holds true on this subsequence (otherwise the 
conclusion is trivial). As we have seen in the proof of statement (2), this implies in 
particular that v e PCd- 

Now let us apply estimate (5.9) to the function We obtain that 

||VGi'=)(..)||^ > ^ ^ ^'^^"^"^ 



4 



1 / JduA^n) 



1-dk + n 1 ^ + [Jdt:,n{'^n)Y 

Jdi+-l,n{Vn) Jdi,n{Vn) 



+E 



^ di+i -di + n 1 V"- ^ + [Jdi+i,n{'"nW n ^ + [JdiA'^n)]'' 

Let us finally pass to the limit as n — )> +oo. Thanks to the jump convergence 
(5.17), the right-hand side tends to || VG^''(f )||2 (assumption (5.23) guarantees that all 
fractions have a finite limit), as defined in Definition 5.6. □ 

5.3 Holder continuity of approximating and limit problems 
Proof of Proposition 2.1 

The existence of a local solution to the system of ordinary differential equations is 
trivial. From now on we identify the vector (ao(^)) ■ ■ ■ ,o,k{t)) with the function u{t), and 
we define the points d^ {i = 0,1, . . . , k, k + 1) as in section 2.4. 

L°° estimate Let t/j E (M) be an even convex function, and let us set 

„1 k 

■^{t) := / ij{u{t)) dx = V((i,+i - dMai{t)). 

Jo .-n 



i=Q 
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Since is nondecreasing, with some computations it turns out that 



1=1 



Let us assume now in addition that il){o-) = if and only if \a\ < HmqUoc- Then we 
have that ^'(O) = 0, ^(t) > as soon as it is defined, and ^' is nonincreasing because 
of (5.29). It follows that '^{t) = as soon as it is defined, which proves (2.15). 

Finite time break- down Let us consider the function 

k 



deD i=l 



The sign of all jump heights is constant as soon as the solution is defined. This 
implies that S{t) is smooth. Computing the time derivative, and rearranging the terms, 
we obtain that (for shortness's sake we drop the dependence on t in the second line) 



S'{t) 



1 1 



di \ai{t) - ao{t)\ 1-4 Wit) - ak-i{t)\ 

k-l 



~^ di+i — di \|aj+i — Ojl {di — Qi-il/ \ l^i — Oj-il — a. 



11 11 

< 



di \ai{t) - ao{t)\ 1-4 \ak{t) - ak-i{t)\ 



< 



1 



\Uo\ 



Since S{t) is clearly nonnegative, this implies that the solution cannot be global, and 
also provides an estimate on the life span. 

Holder continuity up to collision Let us compute the time derivative of the function 
t — >■ G^^{u{t)). Rearranging the terms we obtain that 



i=0 

k 



(t) ai{t) - ai-i{t) 



I f I 1 \ 

with the usual agreement to neglect the two fractions involving terms aiit) with indices 
less than or larger than k. 
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The two factors in each term of the sum are equal due to the system of ordinary 
differential equations. Therefore, the sum can be rewritten both in the form 



1=0 



and in the form 

k 



^ di+i - di \ai+i{t) - ai{t) ai{t) - ai-i{t) 
As a consequence, we have also that 



= ||VGWKt))|| 



d 



Now let us consider the function 



-|G'WKt)) = ||VG'W(^,(t))||^^-^||^'(t)|| 



H{t) ■.= 3kexp(J^G^^\u{t))y 



(5.30) 



(5.31) 



We claim that the descent of H{t) estimates the 1/4- Holder constant of u{t). To this 
end, we begin by computing the time derivative of H{t). Prom (5.30) we have that 

- H'{t) = 2exp (A Gg\u{t))\ \\VG^^\u{t))C WAmT- (5-32) 



Let us estimate the first two terms in the right-hand side. For the first one we have 



that 



l[\ai{t) - ai.i{t)\ 



.1=1 



2/(3fe) 



> min \ai{t) - ai-i{t)f/^. (5.33) 

i=l,...,k 



Let j be the index (or one of the indices) which realizes the minimum. Then from 
Cauchy-Schwarz inequality we have that 



> 



i=0 

E 



1 



ai+i{t) - ai{t) ai{t) - ai-i{t) 
1 



(5.34) 
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Plugging (5.33) and (5.34) into (5.32) we obtain that 

-H'{t) > 2\\u'{t)\&' > \\u'{t)\&' yt e [o,7;i„j. 

Now we integrate in [s, t], and we exploit that H{t) is nonnegative and nonincreasing 
(which follows from (5.32)). We deduce that 

1^ \\u'{T)\&'dT < H{s) - H(t) < H{0) = Skexp (^G^^H^o)) 
for every < s < t < Tsing- Finally, from Holder's inequality we obtain that 

\\u'{r)\\2dr<U ||ii'(T)||f drj \t - s\'/^ 
for every < s < t < Tging- Combining the last two estimates we obtain (2.16). □ 
Proof of Proposition 4.3 

In analogy with (5.31), let us consider the function 



Exploiting equation (2.14) (which is the same as (2.12)), in analogy with (5.32) we 
obtain that 

-Kit) = -2exp(^AGW(^„(i))^(VG'(f)(^„(^)),^;(^)) 

= 2exp [-^G^:\Mt))) \\^G^:Hu„{t))C \KmT. (5.35) 

Let us estimate the first two terms of this product for all t e [0, Tsing,n]- Let Jd,n{t) '■= 
Jd,n{un{t)) denote the discrete jump heights of Un{t), defined according to (2.4). Prom 
(5.5) we have that 

exp (J^G^^\u^{t))^ > exp Qlog (^mm | J,,,(t)| j ^ = (^mm | J,,„(i)|^ ^ . (5.36) 
Moreover, from (5.8) we have that 

||VGi^)(«n(t))||r > [min\JdAt)\j . (5.37) 



Plugging (5.36) and (5.37) into (5.35) we obtain that 

-i/;(t)>2||<(t)||f >||<(t)||f vte[o,r, 



sing,nj 



Now we can conclude by integrating in [s, t] and then applying Holder's inequality 
exactly as in the proof of statement (3) of Proposition 2.1. □ 
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5.4 Well preparation 

In this section we prove our well preparation result (Proposition 4.1). To this end, in a 
time Sn —7- we have three tasks to accomplish: extinguishing all vanishing jump points, 
adjusting the energy, and remaining close enough to Uq. In the next three lemmata we 
examine these three issues separately. Then we make an alternate use of them in order 
to conclude the proof of Proposition 4.1. 

In the following, Jdiuo) denotes the jump heights of Mq, and Jd,nif) '■— Jd,n{un{t)) 
denotes the discrete jump heights of Un{t), defined according to (2.4). 

Lemma 5.9 (Infinitesimal interval convergence) Let D, D' , {u^n}, Uq, Un{t) be 
as in Theorem 3.1. Let {Sn} he any sequence of nonnegative real numbers such that 
Sn ^ as n ^ +00. 
Then we have that 

lim max \\un{t) — UqIU = Q ■ (5.38) 

Moreover, for every n large enough we have that Un{Sn) G PSD"{n),n for some finite 
set D"{n) (which may depend on n) such that D' C D"{n) C D. 



Proof Let w^n G PCn denote the approximation of Uq defined by 

Won (a;) := «o ^ [0, 1], Vn > 1. 

Since uq is already a piecewise constant function, we have that wqu coincides with 
Uq but for the subintervals corresponding to elements of D' . In other words, one can 
think of as obtained by moving every jump of u in the point of the grid on its left. 
In particular, when n is large enough we have that 

Jd,n(«^on) = Jd{uo) Vd e D'. (5.39) 

Of course we have also that w^n uo'm. L'^ as n ^ +oo. 

Let k' \D'\, let Kq be the constant defined in Lemma 5.4, and let us set 

Co := min I Jd(uo) I , Ci := max | Jd(uo)|, C2 ( min < ^ I ) , to: ^^^"^ 



deD' deD' \ ylZ IZ } J Sk'ci 

Let us consider the function 

yn{t) := \\Un{t) - Won||i2((o,i)), 

where the norm is intended with respect to the space variable, and let 

Rn := sup {s > : yn{t) < Ca Vt G [0, s]} . (5.40) 
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Since 

yn{0) = \\uon - WonWl < {\\uon - M0II2 + H^o - WQnh) , 

it is easy to see that y„(0) — > as n — )■ +00, hence i?„ is the supremum of a nonempty 
set when n is large enough. Let us fix uq big enough so that we can apply Lemma 5.4 
for every n > no, and such that 

yn(0) < ^ Vn > no. 

We claim that for every n > no we have that Rn > to, and 

4A;'ci 

yn{t) < yn{0) + -t yt e [0, u]. (5.41) 

Co 

If we prove these claims, then (5.38) is proved. Indeed Sn — > 0, hence Sn < to for n 
large enough, and therefore 

/ Ak'ci \ 

\\Un{t) - Uo\\2 < \\Un{t) - Wonh + ||wOn - M0II2 < ( yn(0) + — Sn j + \\won " Uo\\2 

for every t G [0,5'^]. Since all the terms in the right-hand side tend to zero, estimate 
(5.38) is proved. 

So we are left to prove these claims. Let us consider t e [0, Rn]. Let us apply (5.15) 
with V — Un{t) and w = won (this can be done because the discontinuity sets of u{t) and 
Won are contained in the fixed set D). We obtain that 

min {Ko, \Jd,n{i) ' JdA^on)]} < 6 [y„(t)]'/' < 602/' Vd e D'. 

The right-hand side is less than or equal to Ko/2 due to our definition of C2. Com- 
bining with (5.39) it follows that 

\Jd,nit) - Jdi'^o)\ = \Jd,nit) - JdA'^On)\ < ^C^^^ Vd G D' . 

Therefore, from our definition of co and C2 it follows that 

\Jd,n{'t)\ > \Jd{'^o) \ - \ Jd,n{t) - Jd{uo)\ > Co - QcJ^ > Cq - ^ = ^ (5.42) 

for every d G D' . Now let us compute the time derivative of Unit). Prom equation (2.12) 
and formula (2.8) we obtain that 

y'^{t) = 2y {Un{x,t) -Won{x)) —Un{x,t)dx 

= 2 /' (Unix, t) - Won{x)) ■ uD-^'^ ^ D^^^nix, 

Jo 



l+\Dy^Un{x,t)\' 



dx. 
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Now we apply the discrete version of the integration-by-parts formula (which actually 
is a simple algebraic manipulation of finite sums). We do not have boundary terms 
because D^/"w„(a;, t) is zero both in the last subinterval (1 — 1/n, 1) and in (— 1/n, 0). 
We obtain that 

= 2 \Jd,n{'^On 



deD' 







l + n2| 


l^<i,n(t)P 



The terms of the last sum can be estimated using (5.39), our definition of Ci, and 
(5.42). We obtain that y'^it) < Ak'ci/co for every t e [0,-R„], hence the estimate in 
(5.41) holds true for every t G [0,i?„]. 

Let us assume now that Rn < to for some n > uq. Due to the maximality of Rn we 
obtain that 

/ 7-. \ 4:k'Cl C2 4:k'ci 

C2 = yn{Rn) < yn{0) H Rn < ^ ^0 = C2, 

Co ^ Co 

which is a contradiction. This completes the proof of our claims, hence also of (5.38). 

It remains to understand the location of discrete jump points of ?i„(S'„). First of 
all, the creation of new jump points is forbidden by statement (6) of Theorem B, and 
this proves that D"{n) C D. On the other hand, we have that Un{Sn) — > uq, hence 
statement (1) of Lemma 5.5 imphes that D' C D"{n) when n is large enough. □ 

Lemma 5.10 (Energy adjustment) Let D, D' , {uq^}, Uq, Un(t), he as in Theo- 
rem 3.1, and let k := \D\. Let {Sn} be any sequence such that Sn ^ and 

lim = 0. (5.43) 

n-^+oo 

Let Tsing,„ be defined as in (4-5), and let us assume that 

Tsmg,n > Sn for cvcry n large enough. (5.44) 

Then we have that 

Mm GW(w„(5„)) = G(^)(«o), (5.45) 
where of course G^^{uq) — — oo if D' is strictly contained in D. 
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Proof Since Sn 0, from Lemma 5.9 we know that Un{Sn) Uq. Therefore, the 
Gamma-convergence of Gn^ to imphes that 

nminfG'i^)K(5„))>G«K). 

n— >-+oo 

So we are left to prove the opposite inequality with the limsup, which in turn is 
equivalent to show that 

limsup Gll^\un{Sn)) < M (5.46) 

for every M > G'^^('Uo). 

Let us fix any such M. To begin with, we claim that 

^ El°g (^2 + 4nit)) <M yte [0, S^] (5.47) 

d€D ^ ^ 

for every n large enough. Indeed let us assume that this is not the case. Then there 
exists a sequence {tn}, with tn £ [0, Sn] for every n > 1, such that 

- ^ log - + )>M for infinitely many n's. (5.48) 

deD ^ ^ 

On the other hand, from Lemma 5.9 we deduce that uj^^) ~^ "^q-, hence from the 
jump convergence (5.17) we obtain that Jd,nitn) — > Jd{uo) for every d & D, where of 
course Jd{uo) = if d e D \ D'. In particular we have that 



„ij?<„5E'°«(;^ + ^l»w)=G 

deD ^ ^ 



CO ' 



Mo) < M, 



which contradicts (5.48). 

Prom now on we work in the interval [0, S'^], and we assume n to be large enough 
so that (5.47) and (5.44) hold true. Since Sn < Tsing,n, in this interval we know that 
Un{t) G PSo.n, hence we can use the estimates of Lemma 5.3. Let SQnit) := SQn{un{t)) 
denote the subcritical incremental quotient of Un{t), defined according to (2.5). Let us 
estimate the time derivative of the function t Gn\un{t)). From equation (2.14) and 
estimate (5.7) we have that 

jG^J:\un{t)) = -\\VG(J:\un{t))\\l < -n'[SQn{t)r (5.49) 
for every t G [0, Sn]- On the other hand, from (5.6) and (5.47) we have that 

Gi'\un{t)) < ^l0g(l + [5g„(t)n+^5^l0g(^ + j|jt) 

deD 



71 

< - log {1 + [SQnit)]') +M, (5.50) 
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hence 

[SQr,{t)]' > exp(-[Gl'\ur,{t))-M])-l > -[G('=)K(i))-M]. 
Plugging this estimate into (5.49) we obtain that 

jGi'\un{t)) < -2n [G^:\un{t)) -M] Vt e [0, Sn], 

hence 

G^n\un{t)) < [Gl!:\uon) - M] g-^"* + M G [0, 

Let us estimate Gn\uon)- From (5.50) with i = we have that 

Gl!:\uon) < |log(l + [5Q„(0)]2)+M < ^log2 + M, 

hence 

Gi'\un{Sn)) < ^ • ne-'^'- + M. 

If Sn satisfies (5.43), then the first term in the right-hand side tends to 0. This 
completes the proof of (5.46), hence also the proof of (5.45). □ 

Lemma 5.11 (Discrete jump extinction) LetD, D' , {uQn}, Uq, Un{t) be as in The- 
orem 3.1. Let us assume that D' is strictly contained in D, and let Tsing_„ he the first 
time when a discrete jump disappears, defined according to (4-5). 
Then we have that 

lim Tsing,n = 0. (5.51) 

Proof Let us assume by contradiction that (5.51) is false. This is equivalent to say 
that there exist S > and a subsequence (not relabeled) such that 

Tsmg,n > 5 > for every n large enough. (5.52) 

Let us take any sequence S'„ — ?■ satisfying (5.43), for example Sn = n~^^'^. Due 
to (5.52) this sequence satisfies also (5.44), hence we can apply Lemma 5.10 to this 
sequence. Since D' is strictly contained in D, we have that G^{uo) — — oo, hence 
(5.45) reads as 

lim G^J:\un{Sn)) ^ -<x,. 

Moreover, from (5.44) we have that Un{Sn) € PSD,n for every n large enough. There- 
fore, up to replacing the initial sequence uon uq with the sequence Un{Sn) — > uq (the 
convergence to uq is due to Lemma 5.9), we can always assume that the sequence of 
initial data satisfies (3.1) and 

lim G^^\uon) ^ -oo. (5.53) 

n— >+oo 

32 



Thus from now on we work under this assumption. Let us estimate the time deriva- 
tive of the function t — )■ Gn\un{t)) in the interval [0,Tsing,n]- Using equation (2.14), 
and estimate (5.8), we find that 

Combining with (5.5) we deduce that 

jG^J:\u^{t)) < -exp (-^G(f)(^.„(t))) . 
Integrating this differential inequality in [0, Tsing,n] we obtain that 
I Tsing,n < exp Q G(f)(«n(0))^ - exp Q G(f)(«„(T,i,g,„))^ < exp Q G^fH^on)) ■ 
Thanks to assumption (5.53), this contradicts (5.52), hence it proves (5.51). □ 



Proof of Proposition 4.1 

We argue by induction on k — k', where k := \D\ and k' := \D'\. 

Let us assume that k — k' — 0, namely D — D'. In this case we claim that conclusions 
(4.1) through (4.3) hold true for every sequence S'„ — )■ satisfying (5.43), for example 
Sn '■= n~^^'^. First of all, from Lemma 5.9 we obtain that (4.3) holds true, and Un{Sn) € 
PSD',n when n is large enough (because D = D'). This proves (4.1), and implies that 
assumption (5.44) is satisfied. Therefore, we can apply Lemma 5.10 and deduce (5.45). 
Since k — k', this proves (4.2). 

If /c — A;' > 0, then we begin by applying Lemma 5.11. We obtain the extinction of at 
least one discrete jump in a time Tging.n — > 0. Thus from Lemma 5.9 we have also that 

«n(^smg,n) ^0 (5.54) 

and every element of this new sequence belongs to PSD"{n),n for some D"{n) strictly 
contained in D and such that D"{n) 3 D' . Up to subsequences we can assume that 
D"{n) D" is independent of n. Since the number of possible choices of D" is finite, 
we have only finitely many subsequences to consider. Therefore, it is enough to conclude 
on all such subsequences. 

Setting k" := \D"\, all these subsequences of (5.54) satisfy the same assumptions 
of the initial sequence with D" instead of D, and in particular with k" — k' < k — k'. 
Therefore, the conclusion follows from the inductive assumption. □ 
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5.5 Convergence up to collisions 
Proof of Proposition 4.2 

"Safe intervals" Let Jd{uo) be the jump heights of Wq, and let Jd,n{t) '■— Jd,n{un{t)) be 
the discrete jump heights of Un(t), defined according to (2.4). 

We say that [0,To] (with Tq > 0) is a safe interval if there exists a positive real 
number cq, and a positive integer uq (both may depend on Tq) such that 

\Jd,n{t)\ > Co Vd e D, yt e [0, To], Vn > Uq. (5.55) 

We claim that safe intervals do exist. Indeed let us assume that this is not the case. 
Then there exist a sequence {dk} C. D, a sequence — >■ +oo of positive integers, and 
a sequence — > of positive times such that 

lim Jdu,nu{h) = 0. (5.56) 

re— >'+oo 

Up to subsequences, we can assume that dk does not depend on k. Since tk — > 0, 
Lemma 5.9 implies that Unf,{tk) — >■ uq. At this point the jump convergence (5.17) 
contradicts (5.56). From (5.55) it is also clear that 

^sing.n > Tq for every n large enough. (5.57) 

Boundedness and compactness in safe intervals Let [0,To] be a safe interval according 
to (5.55). In this part of the proof we show that there exist real constants Ci and C2 
such that 

ci < Gl^\unit)) <C2 Vt G [0,To], Vn > 1, (5.58) 

and that there exists v e C° ([0, Tq]; L^) such that (up to subsequences, which we do 
not relabel) 

Un^v in C° ([0, To]; L^) . (5.59) 
Indeed from the monotonicity of the function t Gn\un{t)) we have that 

G^^\ur,{t)) < Gi^H^On) > 0, 

and the right-hand side is bounded from above because of (4.4). Moreover, (5.5) and 
the safe interval assumption (5.55) imply that Gn\un{t)) > fclogco for every n > no. 
This completes the proof of (5.58). Now we exploit a compactness argument. 

• For every t G [0, Tq] (and actually for every t >0) the sequence {un{t)} is relatively 
compact in L^. Indeed from statement (2) of Lemma 5.5 (applied with D' — D) 
we have that 

sup (llltOnlloo + ll-D^^^I^OnlL) < +00, 

n>l 
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hence from statements (3) and (4) of Theorem B we deduce that 

sup {\\Un{t)\\oo + \\D'/''Un{t)\l) < +00. 



n>l 



In other words, we control the L°°-norm and the total variation of Un{t) (as func- 
tions of the space variable). This guarantees the required compactness. 

• The functions ?./,„ : [0,To] — > arc 1/2-Holder continuous, with equi-bounded 
Holder constants, because of (5.58) and statement (2) of Theorem B. Alternatively, 
they are 1/4-Holder continuous, with equi-bounded Holder constants, because of 
(5.57), (4.4), and Proposition 4.3. 

Therefore, Ascoli's Theorem implies that the sequence {un{t)} is relatively compact 
in CO ([0,To];L2). This proves (5.59). 

Passing to the limit in safe intervals Let [0, To] be a safe interval according to (5.55), 
and let v{t) be any limit point of the sequence Un{t). In this part of the proof we show 
that v{t) = u{t) in the safe interval [0,To]. As a consequence, we obtain also that (5.59) 
holds true for the whole sequence, and not only up to subsequences. 

To this end, we write the differential equations in integral form, as in the theory of 
maximal slope curves. Prom statement (1) of Proposition 5.7, we know that equation 
(2.14) implies that (and actually is equivalent to) 

G^:\uM) - Gi'\u^{t)) = \f^ \\u'^{t)\\1 dr+\f^ || VG(^)(«n(r))||^ dr (5.60) 

for every < s < t < To- 

From (5.58) we know that the left-hand side of (5.60) is bounded from above. In 
particular, setting s = and t = Tq, we obtain that 

sup / °i|<(r)||2cir<+oo, (5.61) 

n>l Jo 

sup / ||VG(f)(«n(r))||2 dr < +oo. (5.62) 

n>l Jo 

Prom (5.61) we easily deduce that v e ((0, To); L^), and 

liminf r \K{r)\\ldT > f \\v\T)\\ldT (5.63) 

for every < s < t < To. 

Let us consider now the second term in the right-hand side of (5.60). Due to the 
safe interval assumption (5.55) and (5.57), we can apply Proposition 5.8. Prom (5.24) 
we deduce that 

lim|nf \\VG^:\u^{t))l > \\VG^^\v{t))l W e [0,ro]. 
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Thus from Fatou's Lemma it follows that 




(5.64) 



for every < s < t < Tq. 

Now we consider the left-hand side of (5.60). The functions t — )> Gn\un{t)) are equi- 
bounded and nonincreasing. By the usual compactness result for monotone functions 
(known as Kelly's Lemma, see [1, Lemma 3.3.3] ) there exists a nonincreasing function 
■0 : [0, To] — > R such that (up to subsequences) 

lim Gl^\un{t)) = m Vt e [0,To]. (5.65) 

n— >-+oo 

Now we can take the liminf of both sides of (5.60). Thanks to (5.63), (5.64), and 
(5.65) we obtain that 

> \j\W{r)\\ldT+\f^ pG^S{v{T))g dr. (5.66) 

It remains to characterize the function ^l){t). Coming back to (5.62), and exploiting 
once more Fatou's Lemma, we obtain that 




Therefore there exists a set £^ C [0, Tq] ; with Lebesgue measure equal to 0, such that 
liminf \\VG'^J:\un{t))\l < +oo e [0,To] \ E. 

As a consequence, for every t e [0, To] \ E there exists a (t-dependent) sequence 
Uh ^ +00 such that 

sup||VG(f^)(«„,(t))||2<+oo. 
On this subsequence we can apply statement (2) of Proposition 5.8 and deduce that 
v{t) = lim Un{t) = lim «„, (t) G PCd, 

Jd{v{t)) > Co for every d E D, and 

lim G(f)(«.(t))= lirn G^J^^Mt)) = G^^\v{t)). 

We have thus proved that 

^(t) = G«(i;(t)) yte[0,n]\E. (5.67) 
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The last step is to prove the same equahty for every t G (0,To). To this end we 
remark that v{t) is a continuous function with values in PCd, and is continuous in 
PCd- It follows that the right-hand side of (5.67) is a continuous function. Therefore, 
in (5.67) we have a continuous function and a monotone function which coincide almost 
everywhere in [0,To], hence they coincide everywhere in (0,To). 

Coming back to (5.66), we have proved that 

Gt\v{s)) - G^^\v{t)) > lj\\v'{r)\\ldr + lJ^\\VGg\v{T))\\l dr 

for every < s < t < Tq. Prom statement (2) of Proposition 5.7, this is equivalent to 
say that v{t) coincides with u{t) in [0,To]. 

We have also proved the energy convergence (4.9) for every t e (0,To). 



Continuation up to first jump extinction Let Tqoo be the supremum of all Tq > such 
that [0,To] is a safe interval according to (5.55). Prom (5.57), and the convergence 
results on safe intervals, it is easy to see that 

liminf rsing,„ > Tooo, (5.68) 

lim Un{t) = u{t) e PCd e [0,Tooo), (5.69) 

n— >+oo 

hm G^^\un{t))^G^^{u{t)) Vie(0,rooo). (5.70) 

n— >-+oo 

Let {Rm} be an increasing sequence of positive real numbers such that Rm Tq^o 
as m — >■ +00, and let us set 

Am,n max {\\un{t) - u{t)\\2 ■.0<t< min{i?^, (1 - l/n)T^ing,n}} ■ 
Since [0, R^] is a safe interval for each m, we have that 

Therefore, Lemma 5.1 (standard conclusion) implies the existence of a sequence 
m„ — >■ +00 of positive integers such that Am„,n — >■ as n — >■ +00. We claim that 

Tn :— min |-Rm„, ^1 — — ^ T'sing.n 

is a sequence which satisfies (4.6) through (4.8). Indeed (4.6) is trivial, and (4.8) is 
equivalent to say that ^Imn.n ~^ 0. It remains to prove (4.7). Prom (5.68) we easily 
deduce that r„ — > Tqoo- Therefore, proving (4.7) is equivalent to show that Tqoo = Tling. 
To this end, from (5.69) we immediately deduce that Tqoo < ^ling- Moreover, since 

\\Un{Tn) - m(Tooo)||2 < \\Un{Tn) - m(T„)||2 + \\u{Tn) - m(Tooo)||2, 
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from (4.8) and the continuity of u we have in particular that 



(5.71) 



Let us assume now by contradiction that Tqoo < T'sing, hence m(Tooo) G PCd- In this 
case we claim that there exists 5 > such that [0, Tqoo + 5] is a, safe interval, and this 
contradicts the maximality of Tqoo- 

In order to prove the claim, we argue as in the first paragraph of the proof. If the 
claim is false, then there exist a sequence {dk} C Z}, a sequence Uk — )■ +oo of positive 
integers, and a sequence {tk} of times such that tk G [0,Tooo + 1/^] ioi every k > 1, 
and such that (5.56) holds true. Up to subsequences, we can assume that dk does not 
depend on k, and that tk tends to some limit too £ [0,Tooo]- We can also assume that 
either tk G [0,T„^] for every A; > 1, or tk G [T„j^,Tooc + 1/^] ^'^^ every k > 1. 

In the first case we have that Un^{tk) uitoo) because of (4.8). In the second case 
we have that Un,.(tk) — )■ ^(too) = '?^(^ooo) because of Lemma 5.9 applied to the sequence 
of "initial data" (5.71). In both cases the limit lies in PCd, hence the jump convergence 
(5.17) contradicts (5.56). 

Therefore, we have proved that Tooo — ^sing,n) hence also (4.7). At this point conclu- 
sion (4.9) is exactly (5.70). □ 

5.6 Proof of main results 
Proof of Theorem 3.1 

Global-in-time -convergence We argue by induction on k' — \D'\. li k' = 0, 
then uq is a constant function, and u{t) = uq is the stationary solution. In order to 
prove (3.3) it is therefore enough to show that the function t — )> ||-u„(a;, t) — 'U0II2 is 
nonincreasing. This is true because in this case Unit) — Uq is once again a solution 
of equation (2.12), hence its L^-norm is a nonincreasing function of time because of 
statement (3) of Theorem B. 

Now let us consider the case where k' > 0. In this case we argue in three steps. 

First of all we exploit Proposition 4.1 in order to "well prepare" the sequence of 
initial data. Let {Sn} be the sequence of times provided by Proposition 4.1. We have 
that Un{Sn) — > Wo is a "well prepared" sequence with respect to the discontinuity set 
D', in the sense that the elements of this sequence lie in the corresponding space PSD',n 
when n is large enough, and their A;'-energies converge to the A;'-energy of uq. Now let 
us observe that 



The first term in the right-hand side tends to as n — )■ +00 because of (4.3). The 
second term tends to because 5'„ ^ and u is continuous. It follows that 



max Wunit) — u{t)\\2 < max ||Mn(^) ~ i^olb + max llito — 

tG[0,S„] tG[0,5n] te[0,5„] 




ax \\Un{t) - u{t)\\2 = 



(5.72) 
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The second step is to apply Proposition 4.2 to the "well prepared" sequence of "initial 
data" Un{Sn) — )■ uq. Let {T„} be the sequence of times provided by Proposition 4.2. 
Then (4.8) reads as 

lim max ||m„(<S'„ + — tifOlb = 0. (5.73) 

n^+oo telO,Tn] 

On the other hand we have that 

max ||tin(^) - ii(i)||2 < max WuJSn + t) — u(t)\\2 + max \\u(t) — u(Sn + t)\\2. 

tels„,s„+T„]" te[o,r„] " te[o,T„] " 

The first term in the right-hand side tends to as n — )■ +oo because of (5.73). The 
second term tends to because ^ and u is uniformly continuous. It follows that 

lim max \\uJt) - uitMU ^ . (5.74) 

n^+oo te[Sn,S„+T„]" 

Prom (4.7) we have also that 5'„ + r„ — >■ Tging. Therefore, since 

\\Un {Sn + Tn) - u{Tsing)\\2 < \\Un {Sn + Tn) - u{Sn + T„)||2 + \\u {S^ + T„) - u{Tsing)\\2, 

from (5.74) and the continuity of u it follows that 

Un {Sn + Tn) ^(Tsing). (5.75) 

In the last step we consider this sequence. Due to (4.6), all the elements of this 
sequence belong to PSD',n when n is large enough. On the other hand, the limit 'u(Tsing) 
lies in PCd" for some finite set D" strictly contained in D'. Therefore, we can apply 
the inductive assumption to the sequence of "initial data" (5.75). We obtain that 

lim SUp\\Un{Sn + Tn + t)-u{T,ir.g + t)\\2^0. (5.76) 

Now let us observe that 

sup \\Un{t) - U{t)\\2 < SUp\\Un{Sn + Tn + t) - u{Tsing + t)\\2 

t>S„+Tn t>0 

+ SUp \\u{T^ing + t)- u{Sn + T„ + t)\\2. 
t>0 

The first term in the right-hand side tends to as n — )■ -|-oo because of (5.76). The 
second term tends to because Sn+Tn — > Tging and u is uniformly continuous. It follows 
that 

lim sup \\un{t) -u{t)\\2 = Q. (5.77) 

n-^+oo t>S„+Tn 

Therefore, (3.3) easily follows from (5.72), (5.74), (5.77). 



39 



Global-in-time "uniform" convergence Let Wn{x,t) be defined by 

Wn{x, t) := u {^^, V(a;, t) e [0, 1] x [0, +oo), Vn > 1. 

Since u{t) is already piecewise constant with respect to the space variable, we have 
that Wn and u coincide in /C„. Moreover, both u and Wn are equi-bounded in L°°-norm, 
and for every t > the set Dn{t), where u and w„ could be different, is the union of at 
most k' subintervals of length 1/n. It follows that there exists a constant cq such that 

\\u{t) - Wn{t)\\2 < ^ > 0. (5.78) 

Now for every t >0 we have that 

\\Un{t) - Wn{t)\\2 < \\Un{t) - u{t)\\2 + \\u{t) - Wn{t)\\2- 

Thanks to (3.3) and (5.78), both terms in the right-hand side tend to zero as n — >■ +00 
independently of t, hence 

lim sup \\Un{t) - Wn{t)\\2 = 0. (5.79) 
n-^+oo (>o 

In particular, when n is large enough we have that 

3\\Un{t) - Wn{t)\\T < Ko yt>0, 

where Kq is the constant introduced in Lemma 5.4. Now we exploit once again that 

u and w„ coincide in /C„, and we apply (5.14) with v = Un{t) and w = Wn{t) (their 
discontinuity sets depend on time, but what is important is that they lie inside a fixed 
finite set D). When n in large enough we deduce that 

\\Un{x, t) - U{X, t) ||loo(jc„) = sup \\Un{x, t) - W„(x, t) ||L-((0,l)\D„(t)) 

t>0 

< sup \\Unix, t) - Wnix, t) ||l°°((0,1)) 



t>0 



2/3 



< Ssup \\Un{x,t) - Wn(x,t)\\/2 



t>0 



((0,1))- 



Therefore (3.5) follows from (5.79). □ 



Proof of Theorem 3.2 



Let us consider the double index sequence Am,n '■— Un{m/n) with values in L^. Due to 
(3.6) and (3.7) we have that 

Am,n = Vnim) v{m) 
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Let us apply Lemma 5.1 (refined conclusion) with := V^, and let be a sequence 
such that (5.1) holds true. Let us set T„ := nin/n. Then we have that T„ — )■ and 

For every n large enough we have that the function Un{Tn) lies in some space 
PSD"{n),n, with D' C D"{n) C D. Since the number of these subsets is finite, then 
up to subsequences we can always assume that D"{n) =: D" does not depend on n. At 
this point we can apply Theorem 3.1, with D" instead of D, to the sequence of "initial 
data" Un{Tn) v^o- We obtain that 

lim sup||M„(r„ + t) -'u(t)||2 = 0. (5.80) 
Now let us fix any T > 0. Then for every n large enough we have that Tn <T, hence 

sup \\Un{t) - U{t)\\2 < sup \\Un{t) - u{t)\\2 

t>T t>Tn 

= sup ||m„(T„ + t) - m(T„ + t)||2 

< sup ||'u„(T„ + i) - it(i)||2 + sup ||u(i) - u{Tn^t)\\-2. 

t>0 t>0 

The first term in the right-hand side tends to as n — +oo because of (5.80). The 
second term tends to because u is uniformly continuous. This proves (3.8). 
Finally, (3.9) follows from (3.8) as (3.5) follows from (3.3). □ 
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